Contagion dynamics in time-varying metapopulation networks 
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The metapopulation framework is adopted in a wide array of disciplines to describe systems of 
well separated yet connected subpopulations. The subgroups/patches are often represented as nodes 
in a network whose links represent the migration routes among them. The connections are usually 
considered as static, an approximation that is appropriate for the description of many systems, 
such as cities connected by human mobility, but it is obviously inadequate in those real systems 
where links evolve in time on a faster timescale. In the case of farmed animals, for example, the 
connections between each farm/node vary in time according to the different stages of production. 
Here we address this case by investigating simple contagion processes on temporal metapopulation 
networks. We focus on the SIR process, and we determine the mobility threshold for the onset of 
an epidemic spreading in the framework of activity-driven network models. Remarkably, we find 
profound differences from the case of static networks, determined by the crucial role played by 
the dynamical parameters defining the average number of instantaneously migrating individuals. 
Our results confirm the importance of addressing the time- varying properties of complex networks 
pointed out by the recent literature. 

PACS numbers: 89.75.-k, 87.23.Ge, 05.70.Ln 



Many natural and artificial networks evolve in time, 
with elements that appear, disappear and incessantly re- 
shape their interaction patterns. However, the tempo- 
ral dimension associated to their units and connections 
has started to become accessible only very recently 
thanks to the increasing availability of empirical data [_ 
[Tij . Remarkably, the first analyses of dynamical pro- 
cesses unfolding on these empirical networks have shown 
that the full consideration of the time- varying patterns is 
responsible for a specific phenomenology, clearly distinct 
from the one observed on their time- aggregated counter- 
parts. In particular, the temporal sequence of links and 
their concurrency play a crucial role when the time scale 
of the processes, rn p , is comparable to the time scale of 
the network, tq 0423] ■ This is the case, for example, of 
the spreading of sexual transmitted diseases in a popula- 
tion, ideas in scientific communities, and memes in social 
networks. 

However, the majority of the results obtained over the 
last years hold in two different limits. In the first case the 
network is considered as static. The time scale describing 
the evolution of its structure is much slower than the time 
scale of the dynamical processes, i.e. tq <C tup- In the 
second case instead the networks is modeled as annealed. 
It evolves on a much faster timescale i.e. tq Top, 
and the dynamical processes effectively senses only time- 
averaged properties of the underlying network. While 
these two limits are appropriate for such cases as the 
propagation of failures in technological networks (static 
network) (23j . or to characterize the spreading of easily 
transmitted diseases, such as a novel pandemic (annealed 
networks) [24| . they are not ap prop riate for time- varying 
networks, where t dp ~ r G [3 E HI HI ■ 



The study of dynamical processes in these time varying 
systems, often consists in introducing time windows of a 
certain size At. Thus, the connections between two con- 
secutive intervals are integrated, generating many static 
networks. The dynamical processes are then let evolve 
on top of the sequence of these graphs. This procedure is 
extremely convenient when is possible to define a natural 
time window for the network, but it might introduce un- 
controlled biases due to the size of the considered window 
[2ll ] . For this reason, a different approach that fully ad- 
dresses the time- varying nature of these graphs has been 
recently put forth, and has been adopted to character- 
ize analytically the spreading of infectious diseases and 
random walks [2l|, |22| . It considers activity driven net- 
works dynamics in which the nodes are characterized by 
an activity potential, defining their propensity towards 
creating links. The model reproduces important features 
found in many dynamical systems, being at the same 
time suitable for a simple analytical description of diffu- 
sive processes [2l|, [HI . 

These efforts have considered systems in which nodes 
represent single individuals. Instead, here we address the 
case of time- varying metapopulation networks, in which 
nodes represent units/subpopulations with an internal 
dynamics coupled by the mobility of individuals. In- 
deed, each vertex contains an interacting subpopulation 
connected to its neighbors through links that represent 
routes of migration. Such reaction-diffusion processes 
have been extensively described in static metapopula- 
tion networks f25H32| . but the understanding of the ef- 
fects induced by dynamical connectivity patterns is still 
missing [3(1 [H, Hjl- However, this is crucial to char- 
acterize phenomena as the movements of farmed ani- 
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mals, whose dynamics represents a major issue for pub- 
lic health, and where the evolutionary dynamics of the 
nodes (i.e., slaughterhouses, assembly centers and also 
markets) and links has been shown to be fundamentally 
time dependent [8(. 

We concentrate on contagion phenomena, having in 
mind elementary processes of information or disease 
spreading 35] . For this reason, we consider the prototyp- 
ical Susceptible-Infected-Removed (SIR) model [H, H3] 
for the internal dynamics. The diffusion of individuals is 
modeled by a simple homogeneous diffusion process char- 
acterized by a constant mobility rate p. We consider ex- 
plicitly dynamical connectivity patterns between subpop- 
ulations describing them by activity driven models [2l| . 
In this setting, we describe analytically the early phases 
of the reaction-diffusion process defining the necessary 
conditions for a global spreading of the disease. We ob- 
tain the analytical expression for the threshold value of 
the mobility rate, p* , that determines the outbreak of 
the spreading phenomenon. This result generalizes the 
one concerning static networks, where the same thresh- 
old is determined by the first two moments of the degree 
distribution [29j - [3ll | , and shows that also in the metapop- 
ulation framework the full consideration of the temporal 
patterns of the network is crucial for a satisfactory char- 
acterization of the observed phenomenology. 

The structure of the paper is the following. In Sec- 
tion U we present the basic metapopulation framework. 
In Section |TI] we describe activity driven models that we 
use to characterize the evolution of the links. In Sec- 
tion IIIII within the metapopulation framework we ana- 
lyze reaction-diffusion processes with time varying cou- 
pling patterns. In Section [IV] we discuss the results and 
present the conclusions. 



I. METAPOPULATION FRAMEWORK 

The metapopulation framework is used to describe 
the dynamics of systems characterized by well separated 
but connected units/subpopulations |25l - [33 |. Within this 
paradigm very different systems and processes can be 
described, such as chemical reactions, the evolution of 
populations, and the spreading of infectious diseases. In 
general, any spatially distributed set of subpopulations 
defined by a highly fragmented environment fits this de- 
scription. 

Modeling such systems as networks is natural [29U321 ] . 
Nodes are the basics and isolated units, and links allow 
the movement of particles/individuals. Dynamical pro- 
cesses in metapopulation networks can be described as 
reaction-diffusion processes [3lJ, [32| . While the reaction 
refers to the dynamics taking place inside each node, the 
diffusion concerns the movement of individuals between 
connected nodes (see Figure Q} A for a schematic illustra- 
tion) . In this section we recall the mathematical formal- 
ism needed to tackle these problem for static networks. 



A. Reaction 

As mentioned above, we focus on disease/information 
spreading and contagion processes described using a SIR 
model [3(1, H3| • Individuals are divided in three different 
compartments according to their status, namely Suscep- 
tible (healthy, or unaware), Infected (and spreading), and 
Recovered (cured and no more infective). The evolution 
of the dynamics is defined by two different transitions. 
On the one hand, susceptible individuals in contact with 
infected ones might become infected. This interaction 
happens with a rate f3 and can be described as 

S + I A 21. (I) 

On the other hand, infected people recover from the dis- 
ease spontaneously, with a rate /x described as: 

/Aft (2) 

The inverse of this rate defines the characteristic infec- 
tious period, i.e. the average time each infected individ- 
ual spends in the compartment I. 

The dynamics is determined by the value of these two 
parameters and by the interaction patterns among indi- 
viduals. Here, we consider the simple case of homoge- 
neous mixing. The nodes do not have an internal struc- 
ture, and they all have the same probability of inter- 
action. In particular, whether or not the disease will 
spread locally depends on the basic reproductive number 
Rq = P/ ii, describing the average number of secondary 
cases generated by an infected individual in a fully sus- 
ceptible population [H, HtJ- Diseases characterized by 
Rq < I die out, i.e., the infected individuals are not able 
to sustain the spreading, while for Rq > 1 the disease is in 
general able to take off and spread locally. However, the 
finite size of the subpopulations, and the stochastic na- 
ture of process, may lead to the extinction of the disease 
also in the latter case. In particular, the probability that 
a subpopulation experiences an outbreak give n infected 
initial individuals is 

-Poutbreak 1 ~rm 5 

for all values of R > 1 [H]. 

B. Diffusion 

The diffusion of particles/individuals is a crucial com- 
ponent that drives the entire dynamics of the system. 
Many different diffusion mechanisms have been proposed 
with different levels of realism [H [IJ \M%m- They 
all translate in the definition of a matrix dy characteriz- 
ing the probability that at each time step an individual 
from i travels/diffuses to the subpopulation j. Here we 
consider a very simple diffusion mechanism with the fol- 
lowing properties: 
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• at each time step the probability of traveling ap- 
plies to all individuals in the node independently of 
their origins. These type of processes are Marko- 
vian diffusions. 

• The probability dij is set to be zero if there is not 
a link connecting the two subpopulations. 

• Any individual in a node i, characterized by fcj 
neighbors (degree) will have the same probability 
of reaching any of its neighbors: dij = pAt/ki in 
a time interval At. The diffusion is homogeneous. 
The quantity p defines the mobility rate of individ- 
uals. 

In these settings, we consider a metapopulation network 
characterized by V nodes and N individuals. Ni(t) is 
defined as the number of individuals inside node i at time 
t. We assume that Yli Ni(t) — N for any time t i.e., 
the total population conserved. A convenient description 
of such systems is obtained considering degree classes, 
i.e., by assuming that nodes with the same number of 
neighbors, or degree, are statistically equivalent. 

The average number of individuals in a node of class k 
is then: 



N k 



1 

n 



E 

i\ ki — k 



(4) 



where Vk is the number of nodes of class k. The matrix 
dij can be generalized for degree classes introducing d kk > 
as the probability that an individual diffuses from a node 
of class k to a class k' in an small interval At. Since we 
are considering homogeneous diffusion, we have dkk' — 
pAt/k. 

The variation in the number of individuals in a degree 
class k in a interval At can therefore be written by con- 
side ring th e master equation of a random walk on graph 



d t N k (t) = - P AtN k {t) + kY J N k .{t)d k , k P{k'\k). (5) 



The first term represents the number of people leaving 
the class k. The last term considers the number of 
individuals arriving the class from all other connected 
nodes. In particular for uncorrelated networks, in which 
the probability that a link departing from a node of de- 
gree k connects a node of degree k 1 is independent of the 
destination, the stationary state of Eq. [5] is: 



N k 



(k) 



(N), 



(6) 



where (N) — N/V is the average number of walkers per 
node. The equilibrium describes the final population size 
in each degree class. This results is well know from the 
theory of dynamical processes on networks [35|, 0, |45| . 
The linear dependence on k remarks the effects of the het- 
erogeneous connectivity patterns on diffusion processes. 



C. Invasion threshold 

In metapopulation networks the internal spreading of 
the disease (Rn :> 1), does not necessarily imply a global 
spreading [29l - [33 |. As schematically shown in Figure [T} 
B, in the limit p = the disease is not able to diffuse 
from one subpopulation to the other, but remains con- 
fined to the seeded subpopulations. The spreading at the 
metapopulation level depends on the interplay between 
the time scale of the disease and the time scale of the 
diffusion. The mobility rate p must be large enough to 
ensure the diffusion of infected individuals from seeded 
nodes to others, before the internal dynamics dies out. 
The minimal required value of the diffusion rate p* such 
that for p > p* the disease is able to spread in multi- 
ple nodes affecting a finite fraction of the system, defines 
the invasion threshold in metapopulation networks (see 
Figure [j}B). 

To identify the critical value p* for a general uncorre- 
lated metapopulation network, let us consider a SIR epi- 
demics started from a node of general degree k and size 
Nk characterized by Rq > 1. As discussed in the previ- 
ous section, above threshold the seeded node experiences 
an outbreak with probability given by Eq. [3] In this case 
the number of infected individuals produced during the 
evolution of the epidemics can be written as aN k , where 
a depends on the details of the studied disease. In the 
limit Rq — > 1 for a SIR model, a can be approximated 



by a 



2Cgo-l) 



471 1 . In the SIR evolution, each infected 



individuals stays in the compartment I for an average 
time fj,^ 1 during which it can travel, and potentially in- 
fect other nodes. The diffusion matrix is d kk i , so we can 
write the number of new seeds traveling from a node of 
class k to a node of class k! as: 



A 



OL 

dkk' —N k . 



(7) 



Let us define as the number of diseased subpopulation 
of degree k at generation 0. These nodes are experiencing 
an outbreak at the beginning of the process. They are 
the seeds of the infections. Nodes seeded by one of those 
in the first generation define the set D\ at the following 
generation and so on. If the initial number of diseased 
nodes is small, the progression of the disease can be de- 
scribed using a tree-like approximation relating with 
1 (48l - [50j . Following ref. [3(| we can write: 



D n k 



E 

k' 



-R~ K ' k )P{k\k') 



V k 



This equation considers that each node of degree k' will 
seed the infection in k' — 1 nodes (the total number minus 
the one from which the node got seeded) . The conditional 
probability P(k\k') measures the probability that each 
of the k' neighbors subpopulations has degree k. The 
term 1 — Rq Afe/fc considers the probability to observe an 
outbreak given the number of seeds X k > k arriving from k' 
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to k. The last term in the sum describes the fraction of 
nodes in the degree k still available for seeding. In the 
limit of i?o ~ 1 and in the early stages of the epidemics 
we can simply the equation writing: 



^2 D"- 1 ^' -l)X k 'k(Ro - l)P(k\k'). (9) 



For uncorrelated networks, homogeneous and Markovian 
diffusion dynamics the critical value of p reads: 



(kf 



(fc 2 )-(fc) (R -l)a(N)- 



(10) 



This result clearly shows the strong impact of the topo- 
logical properties of the network on the global dynam- 
ics |29h31| . Indeed, for heavy-tailed networks the ratio 

is typically small, zero in infinite network size, 



implying that the heterogeneities of the metapopulation 
network favors the global spreading of epidemic disease. 
In this case the critical value of p is much smaller than 
the respective value for homogeneous networks. 



II. ACTIVITY DRIVEN MODELS 

In the following, we adopt the activity driven model- 
ing of time varying networks, which represents a natural 
framework to study dynamical processes that evolve at 
the same time scale of the change of the networks [2lj . 
Each node is characterized by an activity potential x, 
describing its propensity to establish links. Formally, x 
is defined as the number of interactions that a node per- 
forms in a characteristic time window of given length At, 
divided by the total number of interactions made by all 
nodes during the same time window. The distribution 
F(x) defines the probability that randomly chosen node 
i has activity potential x, and the latter is bounded in the 
interval e < Xj < 1. In the model, N disconnected nodes 
are initially considered. To each one of them is assigned 
the activity /firing rate Oj = r\Xi, defining the probability 
per unit time to create new contacts or interactions with 
other individuals. The rescaling factor r\ is defined such 
that the average number of active nodes per unit time 
in the system is rj(x)N. In this setting, the generative 
process is defined according to the following rules: 

• At each discrete time step t the network Gt starts 
with N disconnected vertices; 

• With probability a.; At each vertex i becomes ac- 
tive and generates m links that are connected to m 
other randomly selected vertices. Non-active nodes 
can still receive connections from other active ver- 
tices; 




FIG. 1: Schematic representation of a metapopulation system 
and its global threshold. In panel A we show, on the left, the 
diffusion process on the macroscopic level. Red nodes rep- 
resent subpopulations experiencing a local outbreak. Green 
nodes represent subpopulation not yet reached by the disease. 
Nodes with different shades of both colors instead represent 
subpopulation that have been seeded by diseased population 
through the thick links connecting them. On the right we 
show the microscopic diffusion process. Each subpopulation 
contains individuals divided in three compartment according 
to their disease status. Green are susceptible, red infectious 
and black recovered. The arrows represent the flows of mi- 
gration coupling the internal dynamics. In panel B we show 
the critical value of the mobility rate p*. Below this point 
just a infinitesimal fraction of subpopulations is affected by 
the disease. For values equal or larger than the critical point 
the system reach a regime in which a finite fraction of sub- 
populations is affected by the spreading. 



Thus, all interactions have a constant duration tq = At, 
and nodes do not have memory of any previous time 
steps. 

If not specified otherwise, we will consider power-law dis- 
tributions of activity potential i.e. F(x) — Ax 1 that 
reproduce behavior similar to what is observed in real 
data 



21] . At each time step the network is a simple 



• At the next time step t 
network Gt are deleted. 



At, all the edges in the 



random graph with low average connectivity. The accu- 
mulation of connections that we observe by measuring 
activity on increasingly larger time slices T generates a 
skewed Pr{k) degree distribution with a broad variabil- 
ity. Heterogeneities and nodes with a large number of 
connections (hubs) is due to the wide variation of activ- 
ity rates in the system and the associated highly active 
agents. In activity driven models, the creation of hubs 
results from the presence of nodes with high activity rate, 
which arc more willing to repeatedly engage in interac- 
tions. 

The model allows for a simple analytical treatment plj ]. 
The average degree per unit time of a node i character- 
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ized by an activity rate a, will be: 

(fcj)i = ma, + m(a). 



(11) 



The first contribution is due to the m links generated by 
the node when active. The second contribution is due 
to the links that reach node i coming from other active 
nodes. At each time step the average degree per unit 
time in the network will be 



(k)t 



2m(a) 



(12) 



While a snapshot of the instantaneous network would 
show a set of stars, the integrated network, defined as 
the union of all the instantaneous networks at previous 
times, is not sparse. It is possible to show that the degree 
distribution Pr(k) of the integrated network at time T 
takes the form: 



Pr(k) 



1 



Tmrj 



k 



Tmrj 



(13) 



in the limit of small k/N and k/T [2T(. Hence, the inte- 
grated degree distribution has the same function form of 
the activity potential. Fixing an opportune distribution 
for F(x), the model reproduces the connectivity distri- 
bution observed in the data, and considers explicitly the 
dynamical evolution of its structure. 



III. METAPOPULATION MODELS WITH 
DYNAMICAL COUPLING PATTERNS 

In the next sections we present an extension of the 
framework and results obtained in static metapopula- 
tion networks, to systems characterized by time- varying 
topologies. We provide a first modeling approach for 
the study of relevant phenomena as the epidemic spread- 
ing in systems of farmed animals [8], or more in general 
any type of reaction-diffusion process characterized by a 
timescale comparable to the one describing the change in 
the topology of the metapopulation system. 

In particular, we consider V nodes and N individuals 
in a metapopulation network. Each node is characterized 
by an activity potential x extracted from a distribution 
F(x). At each time step, individuals interact locally ac- 
cording to a SIR dynamics and move in or out a con- 
nected node by following a homogeneous and Markovian 
diffusion with rate p. 

Interestingly, due to the effects of time varying cou- 
pling patterns [22j . we find expressions that are com- 
pletely different from the ones observed in the static 
cases. 



A. Diffusion 

The diffusion mechanism has the same properties de- 
scribed in the previous sections, but now occurs on time- 
varying connections. At each time step t, each active 



node creates m random links with randomly selected ver- 
tices, while inactive nodes can be selected by active ones. 
A networks Gt is created and individuals can diffuse fol- 
lowing its links with rate p. A convenient representation 
of such system is obtained considering activity classes. 
Considering the activity rate a = rjx, the number of 
nodes N a in the activity class a at time t are: 



v a 



(14) 



i\ai— a 



where V a is the number of nodes in activity class a. The 
variation of the number of individuals in each class of 
activity in an interval At is given by: 



d t N a (t) = -apAtN a (t) + apmAt{N) 
- pm(a)AtN a (t) + 

+ pAt ( da'a'N a ,(t)F(a') 



(15) 



The first two terms are contributions due to the activity 
of the nodes in class a, while the final two terms repre- 
sent the contribution to inactive nodes due to the activity 
of the nodes in all the other classes. In particular, the 
first term considers that active nodes will let a fraction 
pAt move in others nodes. The second term considers 
that active nodes will connected with m others subpop- 
ulations getting a fraction pAt of their individuals. The 
third term takes into account that active nodes might 
connected with nodes in class a getting a fraction pAt 
of their individuals. The last term considers the number 
of individuals arriving in nodes of class a once these are 
selected by other active nodes. 

As discussed in more details in ref (22j, the stationary 
state of the process is defined by the infinite time limit 
lim t _ s . 00 N a (t) = leading to: 



N„ 



am(N) + (f> x 
a + m(a) 



(16) 



where the constant term 4>i = j daaF(a)N a is the aver- 
age number of individuals potentially moving out from 
active nodes per unit time. In order to better under- 
stand the physical meaning of this quantity let us con- 
sider Eq. [I5l This can be rewritten as: 



d t N a {t) = -F out , a {t) + F in , a (t), 



(17) 



where the two quantities represent the out and in flows 
(negative and positive terms in Eq. I15[) from nodes in the 
activity class a. Since the number of individuals in the 
system is conserved, summing over all activity classes, 
we can write: 

d t J daN a (t)VF(a) = - j da F out , a (t)VF(a) + 

da F in , a {t)VF(a) = 0.(18) 
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The two total contributions must be equal with opposite 
signs. Let us focus on the negative term. The integral 
over all the classes gives the total number of individuals 
moving out of any node in the system: 

F out {t) = J da F out , a (t)VF(a) 

= VpAt J da N a (t)(a + m(a))F(a), (19) 
We can then write: 

F out {t) = VpAt J da N a (t)aF(a)+ (20) 

+ VpAtm J da N a (t)F(a)(a), 

so that: 

F out (t) = VpAtfa + mV P At(a)(W). (21) 

The first term defines the number of individuals moving 
out from active nodes in the system. The quantity <fii is 
defined as average number of people potentially moving 
out of active nodes per unit time. It becomes the actual 
average in the limit p = 1. The second term defines the 
number of individuals moving out from not active nodes 
that have been selected by active ones. 

The stationary state found in Eq. \W\ differs substan- 
tially from the ones found in static ones (Eq. 0, where 
the number of individuals in each node of degree k is 
a linear function of the degree. This is not the case in 
time- varying networks, in fact. The difference is due to 
the effects of the dynamical connectivity patterns that 
play a crucial role in the characterization of the station- 
ary state. In particular, in activity driven networks there 
is a saturation effect for large values of activities. Nodes 
with high activity have on average k ~ m connections 
at each time step, and henceforth a limited capacity of 
collecting individuals (2^ |. 

B. Invasion Threshold 

To better illustrate the results, we focus our attention 
first on the case m = 1, where each active node creates 1 
random link with other subpopulations. We then tackle 
the case of a general value of m. 

1. Case m = 1 

Let us consider a SIR epidemic process started from 
a single node, see Figure [2] for a schematic representa- 
tion. We define D™ as the number of diseased nodes in 
the activity class a and generation n. As seen in the 
previous sections, when the number of diseased subpop- 
ulations is small, it is possible and convenient to adopt a 
tree-like approximation relating Z?" with -D™ -1 . At each 
time step each node might be infected by other seeded 
subpopulations if: 
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FIG. 2: Schematic representation of the macroscopic dynam- 
ics in time- varying metapopulation networks. In the first time 
step, T = 1, just the subpopulation i is diseased (red node). 
It is the source of the infection. By chance it becomes active 
and creates m = 1 random links. Through this connection, 
thick line, some infected individuals to move to the subpopu- 
lation j that will starts to experience the disease at the next 
time step. Next, for T — 2, three others nodes are active. All 
the selected nodes are not diseased so no infected individuals 
diffuse. At the time step, T = 3, node i is active again infect- 
ing another node. The process is repeated for all the other 
time steps, driven by the activity of the each node. 



• The node is active and connects with a seeded sub- 
population that successfully infects it. 

Given X a ' a seeds arriving from a nodes in the class 
a' to the nodes of class a, the destination will ex- 
perience an outbreak with probability 1 — R Xa ' a . 

• The node is not active but is selected by an active 
and seeded subpopulation that seeds it. 

Hence, there are two different contributions: 

D2 = aAtV a J2D:r 1 (l-R Q K '°)(l-^-) v + 

a' \ / 

+ V a J2a>AtD^(l-R^) (l-^U- 

a' \ / 

(22) 

The value X a i a represent the average number of seeds 
arriving from a node of class a' to a node of class a. In 
analogy with Eq. [7] we can write: 

Aa'a = d a , a -N a , ~ pAt-N a ,. (23) 
fi p, 

Indeed, for m = 1 the degree per unit time of nodes that 
can infect or can be infected is 1 + (a) . They have at least 
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one connection, plus a contribution given by the average 
value of activity of all the other nodes. For power-law 
distributions of activity (a) ~ 0(e) -C 1, and we can 
approximate the degree to 1. 

We make further approximations and simplification by 
setting without losing generality At = 1, and by consid- 
ering (i) that for values of Rq close to one (1 — R Aa ' a ) ~ 
A a 'a(-Ro — 1)) an d (ii) that in the early stages of the 
spread the large majority of nodes are not seeded so that 

(1 p — ) ~ 1. We can write: 

D» = aF(a)nY,Dy 1 N a ,+ 

a' 

+ F(a)Qj2 a ' D y 1 ^a', 

a' 

(24) 

where n = p(Rq — 1)^. In general, it is not easy to solve 
directly this equation. In order to find a close form we 
therefore introduce a set quantities: 

a 

c = E ai w 

a 

fa = ^2a h F(a)N a , h=l,2 

a 

(25) 

By multiplying both sides of Eq. (|24| by N a and sum- 
ming, we obtain 

e n = fane' 1 - 1 + niN)^ 1 (26) 

On the other hand, by multiplying both sides of Eq. (pM)) 
by aN a and summing, we obtain 

C = fane 71 - 1 + nfac 1 (27) 

In the continuous time limit, we can write Eqs. (|26p 
and (|27p in a differential form: 

d n e = (fan - l)6 n - 1 + 0(iV)C _1 (28) 

cu = fane 71 - 1 + (fan - 1)^ 1 (29) 

The Jacobian matrix of this set of linear differential 
equations takes the form 

f fan-i n{N) \ 
\ fan fan -1 ) 

and has eigenvalues 

A ( i,a) =0ifi-l ±n^fa(N) (30) 

The epidemic threshold is obtained by requiring the 
largest eigenvalues to be larger than 0, which leads to 
the condition for the global outbreak: 

n(fa + ^fa(N) > 1 (3i) 



we can write the threshold condition on the mobility rate 



(R - l)afa + fa(N) ' 

This equation is the central result of our analysis, and 
it is worth highlighting some aspects of it. First, the 
invasion threshold is not function of the time-integrated 
network representation. It depends on two functions fa 
and fa. As shown in details before, the first quantity is 
the average number of individuals potentially moving out 
from active nodes. The second quantity is an higher or- 
der function of less clear physical interpretation. Second, 
both quantities act together in lowering the value of the 
threshold, meaning that not only the amount of mobility 
observed in the network is responsible for an easier out- 
break of the contagion process, but also its heterogeneity 
act in the same way. Finally, Eq. (|32|) shows that a thor- 
ough consideration of the time- varying nature of the net- 
work is crucial not to introduce biases in the description 
of the processes, also in metapopulation systems. 

We test the analytical prediction of Eq. ([32]) by ex- 
tensive analytical simulations. Considering V = 10 4 , 
m = 1, F(a) = Aa- 2 2 , e = 1CT 3 , (i = 1CT 2 and Rq = 1.1 
we evaluated the number of diseased subpopulation as a 
function of p. Figure |3] confirms that the analytical value 
perfectly matches the simulations. It also shows that at 
large values of p the number of diseased subpopulations 
reaches a maximum, and decreases afterwards. This is 
in an interesting behavior due to the particular choice of 
to = 1 in the activity driven model. When p is large, 
a connected node loses all of its individuals immediately. 
Thus, at each time step very active nodes give to a single, 
and with high probability less active, node all the seeds, 
which then reside for a longer time on the new node. This 
turns out to limit strongly the number of possible new 
seeded subpopulations that can be generated. 



2. General case, m > 1 

The general case of an arbitrary number of links es- 
tablished by active nodes can be derived following the 
same approach described before. The only difference is 
the expression of X a ' a that for active nodes now becomes: 

Xa>a~--K, (33) 
TO /J, 

where we have neglected the contribution to the degree 
given by other active nodes. The value of \ a > a for inactive 
nodes, that are selected by active ones, is the same as in 
the to = 1. Thus, the process is now asymmetric. Active 
and passively selected nodes have in general a different 
degree and a different capacity to accumulate individuals. 
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FIG. 3: The plot shows the number of diseased nodes in a 
metapopulation network, D x , as a function of mobility rate 
p. We consider V — 10 4 populations, (N) = 10 3 individuals 
for each population, m = 1, e = 10 -3 and 77 = 1. We select a 
power-law distribution of activity potential F(x) ~ a; -7 with 
7 = 2.2. For SIR model, we consider /1 = 1CF 2 and Ro = 1.1. 
Each plot is made averaging 2 x 10 2 independent simulations. 
Each one of them started with 1% of random seeds. With 
dots we plot the median of D x . The red triangle represents 
the analytical prediction of the invasion threshold according 
to Eq. [31 




io J 10" 2 ,-. 10"' 10° 



FIG. 4: The plot shows the number of diseased nodes in a 
metapopulation network, D^, as a function of mobility rate 
p. We consider V = 10 4 populations, (TV) = 10 3 individuals 
for each population, m = 3, e = 10 -3 and 77 = 1. We select a 
power-law distribution of activity potential F(x) ~ x ' with 
7 = 2.2. For SIR model, we consider fj, = 10~ 2 and Ro = 1.1. 
Each plot is made averaging 2 x 10 2 independent simulations. 
Each one of them started with 1% of random seeds. Using 
dots we plot the median of Doo . The red triangle represents 
the analytical prediction of the invasion threshold according 
toEq.[3l 



Eq becomes: 

D2 = aF(a)mnY,D , :r 1 N a ,+ 

a' 

+ F(a)flJ2 a ' D u~ lN a>- 

a' 

(34) 

Following the same steps and introducing the same quan- 
tities we get a the threshold: 



(-Ro l)a ( m + + MmNfc + ( m - l) 2 ^ 

(35) 

Remarkably, while this result reduces to the previous 
one, Eq. (|32p . when m = 1, the functional form is now 
different for a general value of m. This is due to the 
possibility of active nodes to release individuals to more 
than one single subpopulation, which interestingly does 
not translate in a simple shift of p*. However, also in 
this case the critical value is a function of the number 
of individuals moving due to active and selected nodes, 
and of its heterogeneity, which again act together in de- 
termining a smaller threshold. 

We test the analytical prediction of Eq. (|35|) by exten- 
sive analytical simulations. Considering V = 10 4 , m = 3, 
F{a) = Aa~ 2 - 2 , e = 1(T 3 , /z = 1(T 2 and R = 1.1 
we evaluated the number of diseased subpopulation as 
a function of p. As show in Figure [4] the analytical value 
perfectly matches the simulations. For values of to > 1 
the behavior observed for large p and m = 1 disappears 
as expected. 



In order to measure more directly the effects induced 
by time-varying couplings, we compare the behavior of 
the fraction of diseased subpopulations, Z?oo, for activ- 
ity driven networks integrated over T time steps with 
the explicit dynamic counterpart we consider so far. In 
the integrated networks links arc accumulated over time, 
and the spreading process is run in the resulting static 
graph i.e. Qi n t = U*Jo^t- ^ s show in Figure when 
the explicit dynamics of the connections among subpop- 
ulation is neglected and integrated over time the critical 
value of the mobility threshold is orders of magnitudes 
smaller. This plot clear show the biases that might be 
introduced by time integrating representations of the sys- 
tem and the importance of considering dynamical cou- 
pling patterns for processes that evolve with comparable 
timescales. These results nicely generalize to metapopu- 
lation systems the findings obtained in networks consti- 
tuted by the interactions of single individuals [2l], [22j . 

IV. CONCLUSIONS 

In this paper we have analyzed reaction-diffusion pro- 
cesses in metapopulation networks characterized by time 
varying coupling patterns, described in the activity 
driven framework. We have focused on the spreading and 
contagion processes of diseases and information, modeled 
through a SIR model. Here, the reaction dynamics in- 
side each node is described considering a homogeneous 
mixing approach, where each pair of individuals has the 
same probability of interact, while the diffusion of indi- 
viduals between connected nodes is Markovian. In this 
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FIG. 5: The plot shows the number of diseased nodes in a 
metapopulation network, Doc , as a function of mobility rate p 
on activity driven network and corresponding integrated net- 
work. The integrated networks are generated by integration 
of activity driven networks over 50 time steps. We consider 
V — 10 4 populations, (iV) = 10 3 individuals for each popu- 
lation, m = 3, e = 10 -3 and r) = 1. We select a power-law 
distribution of activity potential F(x) ~ a; -7 with 7 = 2.2. 
For SIR model, we consider fi = 10~ 2 and Ro = 2.0. Each 
plot is made averaging 2 x 10 2 independent simulations. Each 
one of them started with 1% of random seeds. The dots repre- 
sent activity driven network model and the squares represent 
the integrated model. 



setting, we have characterized analytically the conditions 
for a global spreading as a function of the process param- 
eters and the time varying couplings. Noteworthily, they 
are strikingly different from to the static case, and turn 
out to be a function of the average number of individ- 
uals potentially moving out from active nodes, that in 
the early phases of the process spread the diseases to the 
selected target nodes. These findings not only provide 
the first analytical insight into time varying metapop- 
ulation processes, but clearly point out that dynamical 
couplings affect spreading phenomena also in reaction- 
diffusion processes, and that important biases might be 
introduced by neglecting them. 
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